“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1990-12 


Two nonlinear models of the Free Electron Laser 


Kiel, David H. 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/30622 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


\§ D U DL EY research materials and institutional publications created by the NPS community. 
«iis Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed -- and published -- scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 


NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


-A24 
An inilth 


| 
| 


j rE Ec ~8 aay 
JAN 0 3 1992. 


ye 
a 


THESIS 


TWO NONLINEAR MODELS 


David H. Kiel 


December, 1990 





FREE ELECTRON LASER 


William B. Colson 





Approved for public release; distribution is unlimited 


=0 
=== a 
a= | 
== 
= —_? Thesis Advisor: 
— 
= 
O1 Fy 





-“@ 








Unclassified 


———————— 


SECURITY CLASSIFICATION OF THIS PAGE 


Form Approved 
REPORT DOCUMENTATION PAGE 
UNCLASSIFIED 


2a SECURITY CLASSIFICATION AUTHORITY 3. DISTRIBUTION / AVAILABILITY OF REPORT 



























Vv 


2b DECLASSIFICATION / DOWNGRADING SCHEDULE 


4 PERFORMING ORGANIZATION REPORT NUMBER(S) 


Unlimited Distribution 


5S MONITORING ORGANIZATION REPORT NUMBER(S) 





NAME OF PERFORMING ORGANIZATION 6b OFFICE SYMBOL 


(If applicable) 
33 


7a. NAME OF MONITORING ORGANIZATION 














Naval Postgraduate School 
_ ADDRESS (City, State, and ZIP Code) 


Monterey, Ca 93940,5000 


Naval Postgraduate School 
7b ADDRESS (City, State, and ZIP Code) 






Monterey Ca, 93940,5000 








NAME OF FUNDING / SPONSORING 
ORGANIZATION 


Bb OFFICE SYMBOL 
(if applicable) 


9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 











10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT TASK WORK UNIT 
se ¥ 
TWO NONLINEAR MODELS OF THE FREE ELECTRON_LASER 


ELEMENT NO 
12 PERSONAL AUTHOR(S) ; 
Kiel, David H. 
13a TYPE OF REPORT 13b TIME COVERED 14 DATE OF REPORT (Year, Month, Day) [15 PAGE COUNT 
M : FROM TO De 990 73 
a e De ne 


16 SUPPLEMENTARY NOTATION 








. ADDRESS (City, State, and ZiP Code) 










TITLE (include Security Classification) 

















The views expressed in theis thesis are those of the author and do n 
reflect the official policy or position of the Department of Defense or the U.S. Govt. 


17 COSAT! CODES 
SUBGROUP 


19 A&STRACT (Continue on reverse if necessary and identify by block number) 


t 





18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 











FEL, Macroparticle, Oscillator, 



















The dynamics of the Free Electron Laser are governed by Maxwell's equations 
which causes many highly nonlinear regimes to exist in Free Electron Laser Physics. 
This thesis will examine two such areas and develop simple models to describe 

the highly dynamic and rich behavior two of these regimes. 

In the strong-field, high-current regime, the Free Electron Laser driving 
current can be modeled by a single macroparticle representing the trapped 
electrons. When the trapped electrons act collectively as a macroparticle, 
solutions which include synchrotron oscilations can be found for the self- 
consistent pendulum and wave equations. 

In as FEL Oscillator with low single-pass gain, the evolution of the 
optical wave can lead to sideband development. This phenomenon is studied 
by applying Maxwell's equations to an oscillator with two optical modes and 
20 DISTRIBUTION / AVAILABILITY OF ABSTRACT 21 ABSTRACT SECURITY CLASSIFICATION 
BO UNCLASSIFIEDUNLIMITED =F] SAE AS RPT Unclassified 


22a NAME OF RESPONSIBLE INDIVIDUAL 226 TELEPHONE (include Area Code) | 22¢ OFFICE SYMBO. 
ax 408) 646-2116 p1cw 


DD Form 1473, JUN 86 Previous editions are obsolete SECURITY CLASSIFICATION OF THIS PAGE 
S/N 0102-LF-014-6603 Unclassified 


















O oTic USERS 













i 





SECURITY CLASSIFICATION OF THIS PAGE 


deriving a two-mode wave and pendulum equation. The two-moke wave and 
pendulum equations are implemented numerically on computers so that the 
onset of the sideband can be explored. 





SECURITY CLASSIFICATION OF THIS PAGE 
Unclassified 


OD Form 1473, JUN 86 (Reverse) 











Approved for public release; distribution is unlimited. 
Two Nonlinear Models 
of the 
Free Electron Laser 


by 
David H. Kiel 
Lieutenant, United States Navy 


B.A., University of Colorado 


Submitted in partial fulfillment 
of the requirements for the degree of 


MASTER OF SCIENCE IN PHYSICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1990 


Author: 





David H. Kiel 


Approved by: 








nmdres Larraza, Second Reader 


Department of Physics 


William B. Colsonf Thesis Advisor 


Aecesion For 
NTIS) CRA&! a 


OS Fae 4 








fer betion} | 
feo aie gy Ones 








ABSTRACT 


The dynamics of the Free Electron Laser are governed by Maxwell's equations which causes 
many highly nonlinear regimes to exist in Free Electron Laser Physics. This thesis will examine two 
such areas and develop simple models to describe the highly dynamic and rich behavior two of these 
regimes. : 

In the strong-field, high-current regime, the Free Electron Laser driving current can be 
modeled by a single macroparticle representing the trapped electrons. When the trapped electrons act 
collectively as a macroparticle, solutions which include synchrotron oscillations can be found for the 
self-consistent pendulum and wave equations. 

In an FEL oscillator with low single-pass gain, the evolution of the optical wave can lead to 
sideband development. This phenomenon is studied by applying Maxwell's equations to an oscillator 
with two optical modes and deriving a two-mode wave and pendulum equation. The two-mode wave 


and pendulum equations are implemented numerically on computers so that the onset of the sideband 


can be explored. 
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I INTRODUCTION 


Free Electron Lasers (FELs) are proving to be an exciting and practical source of radiation 
with a wide assortment of applications over a broad range of the RF spectrum. First proposed in 
1970 by John Madey [Ref.1], and experimentally demonstrated by Madey and his colleagues in 1975 
and 1976 [Ref. 2,3}, they are now being considered for use in a wide variety of applications. The 
FEL is based on the original work by Motz who demonstrated that light could be amplified by an 
electron beam propagating through a periodic magnetic field [Ref. 4,5,6]. In 1960, the ubiizon was 
developed by Phillips which operates on similar principals to the FEL [Ref. 7]. 

The FEL has many attributes that make it attractive for a wide assortment of uses. The most 
notable of these attributes is the FEL’s continuous tunability. FELs have been demonstrated to 
operate over a range of frequencies that varied by a factor of 10 [Ref. 8]. Closely related to 
continuous tunability is the "designability" of the FEL. The same basic theory and mechanisms that 
work for a FEL operating in the millimeter range will also work in the tens of nanometer range. 
FELs can also be scaled to high power applications since the gain medium contains only a relativistic 
electron beam in a periodic magnetic field, and does not suffer the nonlinear material affects that 
plague convention atomic lasers in the high power regimes. This also makes FELs very reliable since 
the weakest link in FEL design, the accelerator, is generally considered very reliable. FELs do suffer 
two major drawbacks, size and cost. Some of the current work in FELs is aimed at reducing the size 
which will help with the associated high cost [Ref. 8]. 

Another active area of FEL research is in the high-power regime where they have many 


applications ranging from plasma heating to weaponry. Unfortunately, in the high-power regimes, 











two different areas where nonlinearities arise, and develop simple theories that provide useful insight 
into these dynamic regimes. 

Chapter II will provide the necessary background on classical FEL theory. Chapter II will 
develop a simple theory that treats the bunched electrons in the FEL as a single macroparticle. This 
greatly simplifies the FEL dynamics and allows simple solutions to be derived describing FEL 
behavior in this highly nonlinear regime. The solutions can then be used to study the trapped-particle 
instability, sideband development, and chaos in the FEL. The solutions can also be used for design 
criteria. Chapter IV develops a simple two-mode model of the FEL oscillator by solving Maxwell’s 
equations for two optical modes in the oscillator. The model is then used to study the growth of 
sidebands in an oscillator. A quantitative understanding of this phenomenon will add a tremendous 


degree of control in FEL oscillator design. 


Il FREE ELECTRON LASER THEORY 


A. BASIC FREE ELECTRON LASER PHYSICS 

An FEL consists of two major components, an electron accelerator to produce a relativistic 
electron beam and a wiggler or undulator with a static, periodic magnetic field to “wiggle” the 
electrons with a periodic oscillatory motion. It is this periodic motion of the electrons in the 
undulator that causes radiation in the forward direction. Some of this spontaneous radiation is saved 
in the laser resonator causing stimulated emission that leads to coherent radiation. Since the electron 
beam is highly relativistic, the electrons "see" the undulator’s magnetic field wavelength as being 
much shorter due to relativistic Lorentz contraction. Subsequently, the optical field produced by the 
undulating electrons is of a much shorter wavelength than the undulator’s magnetic field wavelength. 

The above is a Quantum Mechanical description of the basic FEL physics. Despite this, the 
FEL mechanism is essentially classical and can be adequately described with Electro-Magnetic 
Theory. Figure 2-1, illustrates the major features of the classical approach. 

The incoming electrons enter the undulator and begin to oscillate due to the periodic magnetic 
field of the undulator as shown in the top of Figure 2-1. The oscillating action of the electrons causes 
them to emit photons in a manner analogous to a radiating antenna. In the middle of Figure 2-1, one 
undulator wave section is enlarged to show both the optical and magnetic fields of the undulator 
superimposed over the oscillating path of the electron. The forces acting on the electrons are 


governed by the Lorentz force law 


a0B) . _!€ LB x B , y= -“(BeB 2-1 
at mc mc 


where (i = vc and y? = 1 - 8”. The speed of light is c, the mass of the electron is m, the electron 
charge is e, and the velocity of the electron is 1. E and B represent the electric and magnetic fields 


present in the undulator from both the optical wave and the static undulator magnetic field. By 
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ee 2-1: Major components of the Free Electron Laser illustrating the fundamental physics of the 
aser. 





examining the bottom third picture in Figure 2-1, it can be seen that the electrons are effected by the 
optical field and the effect is dependent on the electron’s phase relative to the optical field. This 
coupling between the electrons and the optical wave is best illustrated with an example. When the 
electron has a velocity component both in the z direction and out of the page, and the magnetic field 
of the optical wave is directed downward, the electron will experience a force retarding its forward 
motion and causing the electron to lose energy. Conservation of energy requires that the electron’s 
lost energy go to the optical wave causing subsequent amplification of the optical field. Different 
electron phases with respect to the optical field can also cause the electron to gain energy from the 
optical field. Therefore, careful analysis is required of the electron and optical field coupling to 


produce amplification of the optical field. 


B. ELECTRON DYNAMICS AND THE PENDULUM EQUATION 
In order to fully understand the interaction of the electrons and the optical field inside the 
undulator, the equations governing the evolution of the electrons must be developed. Begin by 


assuming a helically polarized static magnetic undulator field of the form: 

B = B(cos(k,z), sin(kz), 0) , a2 
where B is the field strength, z is the distance along he longitudinal axis of the undulator, andk, 
is the wave number of the undulator’s magnetic field. Also assume that a circularly polarized plane 
wave is present in the undulator cavity of the form 


E, = E(cosy, -siny, 0), B, = E(siny,cosy,0) , 2-3 


where E is the optical field strength and yw = kz - wt + with @ = kc. The wave number for 


the optical field is denoted by k, time by ¢, and the optical phase by 6. 





Substituting equations 2-2 and 2-3 into the Lorentz force equations, 2-1, aad separating the 


electron motion into their transverse and longitudinal components the following is obtained 


dB) 








ae -<[E(1-B)(cosy,-siny,0) +B B(-sink,2,cosky7.0)] , on 
d 
1P) . -£ IER cory - B,siny) + B(B,sinkz - B,cosk,z)] , 
dt mc 
and 
cig = - EIB cosy-B siny] , 2-6 
dt me * y 


where B, = (6, .B , ,0). Since the electrons are highly relativistic, the approximation 1-B, = 0 
can be used to show that E(1 - B.) < BB, [Ref. 8]. This simplifies the transverse equation 2-4 


and allows it to be solved by inspection giving 


i, = - P (coskgesinkyz.0) 2-7 


where K = eBX,/2xmc? is a dimensionless parameter that characterizes the properties of the 
undulator. In a FEL K is usually close to 1. The transverse motion solution, 2-7, assumes perfect 
beam injection so that there are no constants of integration. The transverse motion is on a much 
smaller scale than the average z motion along the z axis [Ref. 8]. 


Substituting the transverse motion, 2-7, into equation 2-6 gives 


q — CAE cos(k,z)cos - sin(k,z)siny)) . 2-8 
mcy 


Using trigonometric identity for the cosine of the sum of two angles and the definition 





€ = (k + k,)z - wt the change of the electron’s energy is found to be 


y= CRE cos(t +o) . 2-9 


ymce 





& is a dimensionless "microscopic" variable that describes the electron’s phase with respect to the 
optical wave and the undulators magnetic field. € is on the "scale" of 2 , the wavelength of the 
optical wave in the undulator, since k>k, and z = ct. The exchange of energy of the electrons is 
easily understood in terms of the electron’s phase from 2-9. For example, electrons with phases such 
that, -7/2 < (¢ + 6) < 2/2 , willreceive energy from the optical wave causing them to accelerate. 
Other electrons at phases, 7/2 $(€ +)<$32/2 will lose energy to the optical wave. Since some 
electrons gain and some electrons lose energy, a beam of initially random, but essentially 
continuously spaced electrons, will begin to bunch during the trip down the undulator. This bunching 
process is the important feature of the FEL gain mechanism as will be discussed in more detail in 
section D. 

To further simplify the dynamics of the electrons in the undulator, it is convenient to relate 


y to € . Using the definition 
yemi-p?=1-B)- Bi , a0 


and recalling 2-7, for the B, motion, we have 


y?=1-f - K’ 2-11 


since Bi=K?/y? . Solving equation 2-11 for B, and differentiating with respect to time gives 
p = UK De 2-12 
Zz y? Y : 
using yY > 1, and B, = . Differentiating the definition of (¢ twice respect to time, and solving 


for B , &ives 





eae 2-13 


Equating 2-11 and 2-12 gives 
tae i 2-14 


eet 


Y Q@+kye i+ 


When one wavelength of light passes over an electron as the electron traverses one undulator 
wavelength, the electron is in resonance with the optical wave, and the electron is oscillating at the 
same rate as the optical field so that € = 0 . Using this, the resonance condition can be derived 
which will be used to simplify 2-14. Starting with the definition of ¢ sii differentiating with respect 


to time yields 
€ =(k+k)cB,-ke =O . 2-15 


Solving 2-11 for B, and substituting into 2-15 yields 


r y? 2-16 





since Y >1. 


Combining equations 2-14 and 2-16 and solving for 7 gives 


es Se eS a 2-17 
(k + kc 2A (kK + ki)e 2k, 


Since kok, , 


ve _ x8 2-18 








Equating equation 2-9, the energy equation, and equation 2-18, and solving for {€ gives 


“ KE 
¢ = oes ——"—cos(f + 6) 2-19 
y?mc 


which is the simple pendulum equation. To put the pendulum equation, 2-19, in dimensionless form, 
define a dimensionless time t = ct/L , such that 0 < t <1, where L is the length of the 
undulator. Applying the chain rule gives d?(/dr? = ((L/c)* so that 


? - 2ew KEL? 4x RELY 


a cos(t + $) = cos(G +o) , 2-20 
y?mc? 


since kc = @, = c2n/X, . N is the number of periods in the magnetic field of the undulator and 
the "°" refers to differentiation with respect to t , as will always be the case. Defining 
Ja| = (4meKELN)/(y?mc?) as the dimensionless optical field strength so that the optical field is 


= |ale’* gives 


g = lalcos(G + o) 2-21 


The pendulum equation given in equation 2-21 is valid for both strong and weak-fields. The 
only restrictions placed on using equation 2-21 are that y be high enough that the electron is not 
bent excessively by the undulator’s magnetic field and that the FEL is not so efficient that y 
decreases substantially as the electron traverse the undulator. The optical field strength |a| in front 
of the cos term explains why stronger optical fields produce rapid electron evolution. 

As the FEL normally saturates, the electrons lose enough energy that they drop out of the gain 
bandwidth. As the electrons continue to lose energy and shift across the gain bandwidth, the beam's 
phase velocity vs é decreases by Av = 4nNAy/y = -2n . By either decreasing the 
undulator’s magnetic field wavelength A, along the z axis, or by decreasing the magnetic field 
strength B along z, or both, the electrons could be accelerated in phase-space to maintain resonance 
with the optical wave. Each method is conceptually equivalent and both are called “tapering” the 
undulator [Ref. 8]. The phase acceleration can be added to the pendulum equation as an added 


constant that is turned on near saturation at time t,. This makes the pendulum equation 





|a| cos(§ +) t<t, 2-22 
5 + |a|cos(f +o) T>T 


where 5 is the necessary added phase acceleration. When the phase acceleration is created by 
decreasing the undulators wavelength, 5 ~ -2nNAA/A, and when the acceleration is provided by 
weakening the magnetic field 5 ~ -4mNK?AB/B(1+K?) [Ref.8]. If 5 > ja] , then the path of 
the electron in phase-space will not repeat itself nor return to its initial conditions. Paths in phase- 
pace that do not repeat themselves are called open orbits. Paths that do close on themselves and 
repeat are called closed orbits. If [a| > 5 some of the electrons can be trapped in closed orbits, 
and with |a| > 6 most of the electrons can be trapped near resonance around € = m producing 


high gain. Successful tapering requires strong optical fields (Ref. 8]. 


C. OPTICAL WAVE EVOLUTION AND THE WAVE EQUATION 

In both FEL oscillators and FEL amplifiers, spontaneous emission rapidly forms a classical 
wave with some degree of coherence since the bandwidth of the spontaneously emitted light is 
comparable to the inverse of the number of undulator periods traversed by the electrons [Ref. 8]. As 
discussed during the pendulum equation’s derivation, there is either energy loss or gain depending 
on the electron’s phase relative to the optical wave. 

The optical wave is taken to vary slowly in time during an optical period (E<wE, $<) 
and vary slowly in space over an optical wavelength (E’<kE, 6’<ko) . This approximation is 
known as the slowly varying amplitude and phase approximation (SVAP) and is consistent with the 
narrow bandwidth expected for any laser. The evolving optical wave can be viewed as a carrier wave 
made up of a single high frequency that is modulated by a complex wave envelope that varies slowly 
in amplitude and phase over many optical wavelengths. It is this complex envelope that meets the 
SVAP approximation. If the complex wave envelope did not meet the SVAP approximation, the light 
from the laser would have a much broader bandwidth than expected and poor coherence. This is 
evident through Fourier analysis. 

To begin the derivation of the FEL wave equation, assume a circularly polarized plane-wave 


10 








is established in the undulator as in 2-3 of the previous section with a vector potential of the form 


~  E,(#) 


A = —4—Giny,cosy,0) , eae 


where all the variables are the same as those defined in 2-3. Using the vector potential in the wave 


equation [Ref. 8] 


where J, is the transverse current from the transverse motion of the electron beam, and the SVAP 


approximation to eliminate all terms with two derivatives, gives 


OE 
. 44 ® <= cos,-siny 0) - 
2-25 


2E = 
AER (sin ,cosyp,0) 7 son, 
c Of c 


The sin and cos functions of w in equation 2-25 rotate fast with respect to the slowly 
evolving complex envelope. To eliminate this fast rotation in equation 2-25, two orthogonal unit 


vectors can be defined, 
¢, = (cosy,-sinw,0) , ¢€, = (siny,cosy,0) , 2-26 


to project into 2-25. This reduces Maxwell’s second-order vector equation, 2-25, into two first-order 


scaler equations 


2). 20 2 
e,{ - 434 a 2 SEG 2-27 
c 


11 








and 


2 2E as 
re (gm iar tes: arenes a 2-28 
c* ar c oO c 

The single-particle current is given by J, = -ecBd°(% - F,) where fF, is the position of the 
i’th electron, and 53(¥ - F) is a three dimensional Dirac delta function [Ref. 9]. This implies that the 
transverse single-particle current with the transverse motion B, from 2-7 of the previous section 
is 
3 e’B . Me ea 
Ja (cos(k,z),sin(k,z),0) 8°(% - F) , 2-29 


ymw, 





where B is the field strength of the undulator’s magnetic field, m is the electron mass and @, is the 
radial frequency of the undulator’s magnetic field. Substituting the single-particle transverse current, 


2-29, into the first-order Maxwell equations, 2-27 and 2-28, and taking the appropriate dot products 








Bives 
_ 2 
J,-e, = eb cos(y + kz8'e- 1) . 2-30 
ymo, 
and 
7 2 
5-8, = “43 siny + kosr@-r) . 2-31 
ymw, 


Summing over all the particles to obtain the total transverse current and averaging both equations at 
a fixed time over a volume element dv that is much smaller than the coherence volume but much 


larger than the optical wavelength, gives [Ref. 8] 
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19E, —_2neK 








= + 2-32 
and 
1 06 ae 2neK 7 ‘ 2-33 
e( 7 pisin(g + >)) , 


where p_ is the particle density and the brackets represent the ensemble average over all the 
electrons. Equations 2-32 and 2-33 can be simplified further by introducing another dimensionless 


parameter 


= 8N(enKL)’p 2-34 
yeme? 
representing the dimensionless current of the electron beam [Ref. 8]. The dimensionless current, j , 


along with the dimensionless field strength lal simplifies equation 2-32 and 2-33 to 


|a| = -j(cos(S + )) , 23) 


and 


= J i + . . 2-36 
$ say 0s )) 


The above two equations describe how the complex optical field evolves with the electron 
beam. With high-current densities, j >, the gain will be high since |a| will be large. Also, 
strong optical fields will slow down the evolution or the optical phase. The above wave equation 
used with the pendulum equation of section B describe the simple FEL optical and electron evolution. 
In general, the pendulum and wave equations have no simple analytical solution, but are easy to 
simulate numerically to show the major features of FEL physics. Section E will present numerical 


simulations that demonstrate important FEL physics. 
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D. FREE ELECTRON LASER GAIN 


1. Low-Current, Weak-Field Gain 
In weak-fields and low-current, an approximate solution to the pendulum equation, 2-21, can 
be obtained by assuming that and |a| are constant. Using conservation of energy in the 
undulator and relating Av_ , from the pendulum equation’s solution, to the change of energy in the 
optical field, an analytic equation can be found to give the gain of the FEL for weak-fields and low- 
current (Ref. 8]. 
The derivation begins by finding a second-order perturbation solution to the pendulum 


equation assuming weak-fields and low-current so that 6 = and |a| = 0 . The solution is 


la, | 
Vo 





C(t) = O,+vot - [cos(C, + vot) - cos(C,) + vtsin(C,)]+... . 2-37 


so that 
la, | 


0 


v(t) = Vv, + [sin(C, +V,t) - sin(,)| + 





a : 1 
S|-zicos + 2v,t) -cos(20,)) + cos(v,T) ~ 


2-38 


1-v,tsin(C,)cos(C, + Vot )] +--- 


At injection, the electrons form an essential continuum that are evenly spread in phase space. 
As the electrons evolve along the orbits in phase-space some will gain energy and move ahead in 
phase space while others will lose energy and fall behind. This bunching occurs on the scale of an 
optical wavelength. Since the optical wave evolves slowly in time and space, each adjacent 2n 
section of phase-space is essentially the same as its neighbors over many optical wave lengths 


allowing equation 2-38 to be spatially averaged over a single 2x length giving 
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(vy) = sf dl, = 


From the resonance condition, it can be shown that Av = 4nNAy/y_ [Ref. 8}. Therefore, 
the change of energy in the optical wave can be found from the average change of energy of the 


electron beam due to conservation of energy. Namely, 


Apts ~(pF dV)ymc7((v) - v4) 2-40 
4nN 


where p_ is the real current density, F is the ratio of the "effective" cross-sectional area of the 
electron beam over the optical beam and is known as the filling factor, and dV is a microscopic 
volume element one optical wavelength long and with the thickness of the optical beam. 

The optical power in a similar microscopic volume with a length of one optical wavelength 
and cross sectional area of the optical beam is P = 2E7dV/8n . Therefore, since gain is defined 


as G = dP/P , the gain is 


[2E7dv/8r] 


eee 


Substituting equation 2-39 for (v) and the dimensionless current density defined in 2-34 into 2-41 
gives 


j{2 - 2 cos(v,t) - v,tsin(v,Tt)] 


G(v,) = ; 2-42 
Vo 


Figure 2-2 is a plot of the low-gain curve. It is interesting to note that the optimal gain 
occurs at V, = 2.6 rather than on-resonance or far off-resonance. This is because near resonance, 
the electrons can not drop a large distance through phase-space. Far above resonance, Vv, > 1 the 
gain will oscillate many times along the undulator resulting in small final gain. In this low-current, 


weak-field regime, gain is not exponential. 
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Figure 2-2: The weak field gain spectrum for low current. 


2. High-Current, Strong-Field Gain 

A gain spectrum, G(q,,v,) , is easily obtained numerically by calculating the gain at the 
end of a simulation at each point (V,,a,) with current j. Figure 2-3 is the result of such a 

calculation. 
For weak-fields |a,| S m_ the profile in Figure 2-3 is essentially that derived analytically 
in section 1 and plotted in Figure 2-2; ie it is anti-symmetric with a peak gain of G = 0.135; at 
Vv, = 2.6 . However, as the initial field strength increases, the peak gain of Figure 2-3 shifts to 
higher v, . This is due to the large initial separatrix, whose height is 2V\a| , that can now trap 
more particles farther off-resonance. The trapped electrons can then lose more energy dropping 


through phase-space. Also, the peak gain diminishes with higher field strengths. This is saturation, 
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Figure 2-3. The gain spectrum for high current with increasing initial field strength. 


and is caused by the trapped particles in the closed orbits of phase-space. After trapping, the particles 
“slosh” back and forth in phase-space producing synchrotron oscillations as will be described in detail 
later in section E.4. Finally, note that the gain spectrum distorts and the gain bandwidth broadens 
in v, [Ref. 8]. 
3. Gain Degradation Due to Beam Quality 

Free Electron Laser design often involves a trade off between high-current and beam quality 
from the accelerator. Since the optical wavelength is usually set by the application, the dimensionless 
current j is determined by either the real current or the number of undulator periods. Increasing the 
real current degrades beam quality by producing a beam with more energy spread which degrades the 


maximum theoretical gain. Conversely, increasing the number of undulator periods to increase 
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coupling, may decrease the maximum theoretical gain since the undulator is now more sensitive to 


angular spread of the beam. 
Emittance is the measure of the electron beam quality, and is a constant given by the product 
of the average radial and angular spread of the electrons €, = yo, . It can be shown (Ref. 8] that 


the average decrease of phase velocity due to the emittance of the beam is 


Av, ae 2nN 
1+K? 





(K*koyo + y70;) . 2-43 
In order to get an understanding how the spread of phase-velocities affects gain, it is 
necessary to assume weak-fields and combine the wave and pendulum equations. First, solve the 
pendulum equation for the i’th electron with a first-order perturbation €, = € + vt + 0 , and 
substitute this solution into the wave equation. Integrating over all initial phase velocities, fav, : 


gives an integro-differential equation governing the evolution of the optical field 


ate) we. f ‘dt! F (te "a(t -t) 2-44 
2 vo 


where F(t’) = fea fig)e “is the characteristic function for the distribution fig) , and fig) 
is the initial distribution of electron phase velocities v, = v,+qg about v, . In the integral 
equation, 2-44, all reference to electron phases has been removed, but there still exists the average 
distribution of initial electron phase velocities (Ref. 8]. 

Two common distributions exist, the Gaussian 


2 
f,(9) = iia La 2-45 


y2n 6, 


and the exponential 


exp(q/o,) 


S,(9) = for q<O . 2-46 


18 





These common distributions have an easy functional form, but real beams generally have additional 
factors complicating their shape and their exact functional form is unknown. Still the Gaussian and 
exponential forms provide realistic insight to the gain degradation problem which can be explored by 


substituting one of their functional forms into the integral equation 2-44 and integrating it numerically 


1n (1+G) 


AE 
VA tA vA, 


g 





Figure 2-4: The weak field gain profile for low current with an increasing Gaussian energy spread. 


to obtain the gain profile in Figure 2-4. Several features are immediately evident from Figure 2-4. 
First, the peak gain decreases steadily as the beam quality decreases with bigger o . The peak gain 
also shifts further off-resonance with poorer beam quality and finally, the gain bandwidth becomes 
wider. These general features are common to many other beam distributions, but the details of the 
degradation can be dramatically different. With high-currents, j>1 , and poor beam quality, the 


gain is exceeding difficult to predict and can be off by several orders of magnitude (Ref. 8]. 


19 











4. High-Current, Weak-Field Gain 

In the high-current, weak-field regime, the integral equation, 2-44, can be used to develop 
analytic solutions for the optical field strength and phase. In order to solve the integral equation, 
consider an FEL with perfect beam quality, f(t’) = 1 , and starting on-resonance, v, = 0 . 
Taking successive derivatives of 2-44 then gives, a=i ja(t)/2 , which has a general solution 
inv-lving three roots of the form a,exp(a,t) . The coefficients are determined by the initial 
conditions a(O) = a, and a(0) = a(0) = Q . Theexponents a, are complex roots of the cubic 


equation a; -ij/2 = 0 giving a solution which is the sum of all three roots 


a(t) = Flexp2)'"G + (B)t/2 + exp(i2)"i-V3)t2 + expt]. 2-47 


During the bunching process, t <1, = (2/j)'? < 1 , 2-47 reduces to 
a(t) = a1 + ijt?/12 +...) , 2-48 


implying that the optical field strength remains esser...a.y constant while the phase increases 
proportional to 1? . With high-currents only the fastest growing root dominates after the bunch 


is complete giving 


la(t)I = exp i243 t/2, lt) = (i2)""x/2., 2-49 


implying that the gain is exponential until the onset of strong-fields. With large currents, j>1 , 


the optical phase shift is large as well. 


FE. NUMERICAL SIMULATION OF FREE ELECTRON LASER DYNAMICS 
The pendulum and wave equations, 2-22, 2-35 and 2-36 are three nonlinear differential 
equations that have no easy, general analytic solution. However, they are easy to implement 


numerically and can be used to illustrate major features of FEL physics. 


1. Electron Bunching 

Figure 2-5, illustrates electron bunching with a monoenergetic, low-current beam with a 
weak optical field present. The electrons are at resonance, vo = g = 0 . The large plot in Figure 
2-5 shows the evolution of the electrons in phase-space. The phase-space plot is a plot of the 
electron’s phase, €, versus the electron’s phase-velocity and reveals many of the FEL’s dynamics. 
The initial positions of the electrons are shown in light grey and become darker as they evolve. The 
separatrix separates orbits that repeat and close on themselves, closed-orbits, from open-orbits that 
do not close on themselves. The separatrix shifts position with time due to the term in the 
pendulum equation. The separatrix is plotted in this simulation as the black line in phase-space. It 


can be shown [Ref. 8] that the separatrix height is 2 ja|'? . The remaining two smailer plots show 
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Figure 2-5: Electron bunching with sample electrons on resonance. 
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gain and the evolution of the optical wave’s phase @ . Attime tT =" the electrons enter the 


undulator evenly spaced in phase with respect to the optical wave, and begin travelling down the 
undulator slightly slower than the optical wave. The electrons whose phase m/2 < © + @ S 3n/2 
lose energy and begin to fall back in phase space while those electrons with -n/2 < € + o <n/2 
gain energy and accelerate in phase space. This bunching, under the right conditions, can trap 
particles in the separatrix and leads to substantial gain and subsequent saturation as will be 
demonstrated in section E.3. In this example, the gain is minimal since approximately half the 
electrons gained and half lost energy. The little gain that exists is due to the optical phase shift,  , 
that causes the separatrix to shift slightly from its initial position leaving the average position of the 
electrons slightly greater than 1/2 . 
2. Low-Currents 

To achieve non-trivial gain, it is best to start the electrons off-resonance as illustrated Figure 
2-6. Here the electrons were started at v = 2.6 and a, = 2 with alow-current and weak initial 
fields. More electrons were used to show the virtual continuum of electrons present in the undulator 
as well as increase the accuracy of the simulation. Note that the most dense bunching occurs near § = 7 
where the phase space paths are directed downward implying the highest loss of energy for the 
electrons, and henceforth the largest gain. Also, note that the electrons followed the open orbits 
outside the separatrix, and that there was only a minimal shift in the optical phase. 

3. Strong-Fields, Low-Currents, and Saturation 

Figure 2-7 shows a simulation that demonstrates saturation in a low-current, but strong-field 
FEL as might occur during a single pass in a FEL oscillator. In this simulation the electron beam 
is no longer monoenergetic, but has a Gaussian distribution of phase-velocities with a width of 
0, = 5 . This more realistic situation prevents all the particles from being trapped inside 
separatrix, which ultimately lowers the maximum theoretical gain that would have been obtained with 
a monoenergetic beam as discussed in section D.1 and D.3. The dots in phase space represent the 


final positions of the electrons at the end of the simulation rather than their time evolution through 
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Figure 2-6: Starting electrons off resonance increases gain for a single pass through the undulator. 


phase-space. The electron dynamics of this simulation illustrate several important FEL features. Due 
to the initial energy spread in the beam, many of the electrons did not become trapped inside the 
separatrix and are more spread out in phase-space than those that are bunched inside the separatrix. 
Inside the separatrix, the electrons bunch quickly due to the strong-fields and the optical power grows 
exponentially as the electron bunch moves past (+) = m at about t = 1/2 . It is the peak to 
peak height of the separatrix, 4 |[a|’” , at this point that determines the amount of energy that the 
electrons can transfer to the wave since Av(t) = 4 |a|'® . This change of phase velocity is enough 
to move the electron bunch across the gain bandwidth. The growth begins to slow down until near C+ = 1/2 

where the electrons lost all the energy possible, and flow into the -2/2 <Q + 6 <n/2_ region of 


phase space where the electrons will take energy from the optical wave anu cause the drop in gain. 
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Figure 2-7: A low current beam in strong fields at saturation. 


4. High-Current, Strong-Fields and Synchrotron Oscillations 
Figure 2-8 corresponds to a high-current beam with j = 10000, the optical wave grows 
exponentially, as discussed in section D.4, until the onset of saturation. The saturation occurs with 
similar dynamics as the example of section E.3. However, after saturation and the downward trend 
of the gain evolution caused by the electron bunch flowing into the -2/2 < (+ < n/2_ region 
of phase space, the bunch continues to be driven by the optical wave around the closed orbit inside 
the separatrix and back into the m/2 < (+6 < 3n/2 region of phase space. There the electrons 
transfer their energy back to the optical wave which causes the second increase of the optical power 
that is seen in the gain evolution. This oscillation in the gain evolution is known as a synchrotron 

oscillation. Also note how rapidly the optical phase evolves with time. 
In Chapter I, it will be shown that the frequency of the synchrotron oscillation is 
jin 


v, = la . Chapter HI will also develop a simple theory using the concept of the a single 
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Figure 2-8: High currents and strong fields cause synchrotron oscillations due to the electron bunch 
“sloshing” around phase space following the closed obits in the Separatrix. 


“macroparticle" to represent the electron bunch. This will allow the development of an approximate 
solution of the wave and pendulum equations in this highly nonlinear regime. 
5. The Tapered Undulator 

This example illustrates the effect of tapering an undulator to achieve higher gain. At 
injection, the electron beam had a Gaussian energy spread with a width of 6, = 5 , anda current 
of j = 2000. The tapering begins near saturation at Tt, = 0.6 with a tapering rate of § = 702. 
Note that most of the electrons were trapped inside the separatrix, but that some remain outside the 
separatrix and are randomized in phase-space. These randomized electrons play a minimal role in 
the gain mechanism whereas the trapped electrons bunch near the critical point lose energy, and lead 
to gain. The "sloshing" of the electrons due to synchrotron oscillations is still visible in the gain 


evolution, but now the closed-orbit area of phase-space is relatively small due to the tapering of the 
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Figure 2-9: A strongly tapered, high-current undulator. 


pendulum equation. This causes the electrons to bunch near (+ = 3n/2 so that there is 
continual growth. The electrons still “slosh”, but over a much smaller area of phase-space which is 
why synchrotron oscillations play a relatively small role in strongly tapered undulators compared to 


untapered undulators. 
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il MACROPARTICLE MODEL 


A. INTRODUCTION 

The high-current, strong-field regime of the FEL is highly nonlinear and difficult to analysis 
quantitatively. However, simple analytic models can be developed by treating the particles trapped 
in the closed orbit region of phase-space as a single macroparticle. The macroparticle approach was 
originally developed and used to derive simpler expressions of the pendulum and wave equations that 
could be used to facilitate FEL design (Ref. 10,11]. The macroparticle approach has also been 
successfully used to study the stability of Klystron-type FEL’s [Ref. 12], and to explore the limits of 
the pendulum model for FEL dynamics [Ref. 13]. In this chapter, simple solutions to the wave and 
pendulum equations are derived using the macroparticle approach. The solutions allow quantitative 
analysis of the optical wave evolution, and provide useful insight into the trapped particle instability. 

The macroparticle model presented here is first solved for the tapered undulator where 
synchrotron oscillations are small. Then a slightly more complex model is developed leading to 
solutions that describe synchrotron oscillations in the untapered undulator. Comparisons of the simple 


theory with particle simulations will be interspersed to demonstrate the validity of the model. 


B. MACROPARTICLE MODEL 

In the high-gain, strong-field regime, the trapped electrons bunch in the closed-orbit region of 
phase-space allowing them to be treated as a single macroparticle. The untrapped electrons will 
randomly disperse in phase-space and play a greatly diminished role in the gain mechanism. This 
removes the required averages from the wave equations and simplifies the system, equations 2-22, 


2-35 and 2-36, from 2n + 1 equations to a system of 3 equations, 


C= 8er -%,) + la[cos(C +o) , 3-la 
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lal = -j,cos(6 +) , 3-1b 
lalo= j,sinS +o) , 3-1 


where Q(t - t,) is the step function for the tapering constant, 5. The trapped particles are described 
by a reduced particle current j, = 0.6; . The macroparticle remains trapped near resonance, at the 


bottom of the potential well, so that t= 0. 


1. Solutions Without Synchrotron Oscillations _ 
The near-resonance assumption allows for the immediate solution of the pendulum 
equation, 3-1la, to give cos( +) = -5/|a| . When combined with the wave equations, 3-1b and 


3-1c the result is 


pile? om 
Ja] 
R 3-2b 
: §? : 
lal¢ = isto = Jr 


A properly designed undulator requires §?< lal? to insure that there will be closed-orbits that trap 
the electrons which is consistent with the approximation in 3-2b. The solution to 3-2a is found to 


be 
Ja(t)|? = 2j,8(t-t,) + Ja,/? . 3-3 


The solution to 3-2b is found by combining 3-3 with 3-2b and integrating with the appropriate initial 


conditions yielding 
(t) = Ja, |[(1 +2j,5(t-t,)/ la, [?)'? - 11/8 3-4 


Without tapering the solutions, 3-4 and 3-5, simplify to 





la(t)|? = la? , O(t) = o,+ 2 t. 3-5 


a, |? 
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Figure 3-1 illustrates a simulation of a tapered undulator with an initial optical field 

a, = 10 , a dimensionless current j = 20000 , and a tapering rate of 5 = 270m . The 
electrons are initially placed with a random Gaussian spread in phase velocity related to an energy 
spread in the electron beam by 4% NAy/y = Av. The random Gaussian distribution about v, is given 
by f,(@) = exp[-g?/207, Wy2n O,, where the i’th electron’s phase velocity is v, = v, +q. In Figure 
3-1 the width of the distribution is 6, = 3 . The large left window in the figure is a plot the 
electron phase-space; electron phase, ©, versus electron phase velocity, v, where each dot is the 
electron’s final position. Note that the trapped electrons are bunched inside the separatrix while the 
untrapped electrons were dispersed randomly. The separatrix separates the closed-orbit region of 
phase-space from the open-orbit region and is plotted as the thin black line in phase-space. The 
electrons inside the separatrix are bunched near x which causes the continual growth of the power. 
Synchrotron oscillations are visible as small modulations in the gain evolution and are caused by the 


"sloshing" of the electrons inside the separatrix. The top-right window plots the power of the FEL 
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Figure 3-1: High Current, strongly tapered undulator. 
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over time with the derived analytic solutions superimposed in grey. The simulation’s gain evolution 
is plotted in black. The flat initial portion of the analytical solution, t = 0 , and the exponential 
gain portion, t, < t < T, , are from the results stated in Chapter I, section D.4. The final portion 
of the analytical solution is the macroparticle solution that begins after saturation when tapering 
commences at t = 0.3. Sixty percent of the electrons were assumed trapped in this example. After 
saturation when tapering commences, the power grows linearly. The bottom-right window plots the 
evolution of the optical phase with the same three regions of the analytical solution superimposed in 
grey. Note that the phase grows linearly in time during the exponential growth as predicated. There 
is close agreement between the analytic solution and the simulation. A gainof G = 2x10 is 
achieved. 

Figure 3-2 illustrates a tapered undulator with a lower current of j = 10000 anda 
smaller taper rate of 65 = 200n , but a strong initial optical field a, = 100 . The initial electron 


beam has a Gaussian spread i: pnase velocity with a width of 6, = 5 . Saturation occurs at 
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Figure 3-2: High initial field in a strongly tapered undulator. 


30 





t, = 0.2 since the initial field is near saturation. The gain grows linearly after saturation. Figure 
3-3 is similar to the simulation of Figure 3-2, but the initial field is much weaker, a, = 10 . This 
delays saturation until t, = 0.4 . Both simulations assume j, = 0.6j. In both figures, the simple 
macroparticle model compares extremely well with the simulations. Synchrotron oscillations are 


visible in both figures as the small modulation in the gain and phase evolution. The analytic solution 
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Figure 3-3: Weak initial field in strongly tapered undulator. 


can not follow the synchrotron oscillations in the gain evolution since the analytic solution assumed 
that the macroparticle was rigidly trapped at the bottom of the optical potential well. However, the 
macroparticle can oscillate slightly as is evident by the closed orbits inside the separatrix and the 
slight modulation of the gain evolution. With tapered undulators, these oscillations are relatively 
small and can be ignored. In untapered undulators, synchrotron oscillations are substantially more 
pronounced and must be accounted for. The next section will account for the oscillations of the 


macroparticle in the potential well and accurately predict the synchrotron oscillations. 
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Figure 3-4 shows a simulation with parameters similar to the LLNL ELF experiment, but 


with better beam quality (Ref. 14]. The dimensionless current is j = 7000 with a random phase- 
velocity spread of width o, = 3 . The initial optical field is a, = 12 , and the taper rate is 
5 = 80x . Again, the analytic solution closely follows the simulation. The simulation and the 


analytic solution both show efficiency of about 30% which compares well with the experiment. 
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Figure 3-4: A simulation of an FEL with parameters close to LLNL ELF is compared with the 
macroparticle model. 


2. Solutions with Synchrotron Oscillations 
The macroparticle model can be extended to include synchrotron oscillations for untapered 
undulators, 56 = 0, where the effect of synchrotron oscillations is more pronounced. With the 
assumption that the macroparticle is trapped at the bottom of the optical potential well, ‘a = 0 , the 


pendulum equation, 3-la, can be solved for an equilibrium position 
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2. tt 3-6 
Cat 9 
The macroparticle is now allowed to oscillate near the bottom of the optical potential well, around (° 


so that 


C= Great) = 5 -o +x), 3-7 


where x(t) < m™ . Assuming that the optical field does not change rapidly after saturation, 
{a| = $ = 0 , the solution for x(t) can be found by combining 3-7 with the pendulum equation 


3-la. Using a small angle approximation for x(t) yields 


@ =X = |a,|cos(G +x(t) +) = -la,[x(t) , 3-8 


where |a,| was substituted for |a| since |a| is nearly constant after saturation. Equation 3-8 has 


a solution of 
x(t) = Ccos(v,(t -t,)) . 3-9 


The frequency of the synchrotron oscillation is v, = |a,|'* . The initial conditions for the 
macroparticle in phase-space aie determined phenomenologically by comparing with simulations. The 
best solution is found to be €, = 0.6 and v, = 0 as expressed in 3-9. 


Combining 3-7 with 3-1b yields 


Ja|= ~j, cos(* + x(t) +9) 5) 


= jpx(t) , 


after the small angle approximation is applied. Combining 3-9 and 3-10 produces a solution for the 
optical field strength 
= Ink 


s 


sin(v,(t -7,)) + la,|] . 3-11 





33 





Combining the optical phase expression of 3-2b with the new optical field strength solution from 3-11 


gives 
dr ? e ' 
o(t) = >, + Te ~T,) + ;cos(v,(t-%,)) , 3-12 
s 


ja,|?v 





sz 


after expanding to first-order in x(t) and integrating. 

Figures 3-5 and 3-6 compare the improved macroparticle model for untapered undulators 
with simulations using current j = 1000 , and initial optical fields of |a,| = 1 and |a,{ = 30 
respectively. These examples illustrate the macroparticle model’s excellent ability to predict the onset 
of saturation. Both simulations use nearly monoenergetic beams described by 06, = 1 with 1000 
sample electrons and 300 time steps. Note that in Figure 3-5, the FEL just saturates, whereas in 
Figure 3-6, the FEL saturates and undergoes nearly two synchrotron oscillations. This is due to the 


high initial field that causes saturation to happen at a much earlier time. 
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Figure 3-5: A comparison of the macroparticle model's solutions with a weak initial field simulation 
and current j = 1000. 
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Figure 3-6: A comparison with higher initial fields and a current of j = 1000. 


Figures 3-7 and 3-8 demonstrate the macroparticle model capabilities with much higher 
currents of j = 10000, and initial fields of {a,| = 10 and |a,| = 100 . Both simulations have 
nearly monoenergetic beams, 0, = 1, with 1000 sample particles, use 600 time steps, and assume 
sixty percent of the electrons are trapped. In both examples, the macroparticle model predicts 
saturation very well and follows the synchrotron oscillation. The analytic optical phase solutions also 
follow the simulation closely and show the synchrotron oscillations. With the higher initial field, 
saturation occurs sooner. 

Figure 3-9 shows a simulation of an untapered undulator with a much higher current of 
J = 20000, and nearly perfect beam quality, 6, = 1. Note the very tight bunch inside the separatrix 
even after nearly four synchrotron oscillations. Once again the analytic solution follows the 
simulation closely, but the predicted synchrotron frequency is slightly higher than in the actual 


simulation. "ote that the analytic power solution, P(t) , after saturation is sinusoidal but in the 
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Figure 3-7: High current with 3 synchrotron oscillations. 


simulation the synchrotron oscillations are not exactly sinusoidal. This is due to the assumption that |a| = 0 
after saturation which is not completely accurate. In the simulations, the peaks of the power curve 
are sharper than the troughs since |a(t)| is not completely constant as was assumed. In Figure 
3-10, Coulomb affects have been examined for this example and appear to have little effect on the 
theory. . 

As a final example, Figure 3-11 shows the comparison between the solutions of 3-7 and 
3-8 and a simulation using 5000 electrons that has parameters similar to the LLNL ELF untapered 
undulator {Ref. 15]. The initial field strength is la, | = 10 , the dimensionless current is 


j = 3350 and the random phase-velocity spread is o, = 3 . From the power evolution, P(t) , 
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Figure 3-8: Higher initial field causes saturation to occur sooner. 


it is apparent that the electrons undergo about two synchrotron oscillations, and that the analytic 
solution closely follows the simulation. The analytic phase solution also shows synchrotron 


oscillations as it follows the simulation. 


C. CONCLUSION 

The macroparticle model provides simple analytic solutions to FEL dynamics in the highly non- 
linear regime with high-currents and strong-fields. The results could provide the necessary tools to 
lead to multimode theoretical studies of the trapped-particle instability, sideband development and 
chaos in the FEL. Such simple analytic results can be used in the design of future high-power, high- 


gain FEL’s. 
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Figure 3-9: A high current example, j = 20000, with four synchrotron oscillations. 
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Figure 3-10: A high current FEL with many synchrotron oscillations compared to the macroparticle 
model with Coulomb affects considered. 
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Figure 3-11: Macroparticle model compared with untapered undulator that has parameters close to 
the LLNL ELF FEL. 
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IV TWO MODE THEORY 


A. INTRODUCTION | 

In the FEL oscillator, the optical evolution can develop sidebands and become highly chaotic 
over many low gain passes just as the FEL amplifier can develop sidebands and chaotic behavior over 
one high gain pass. This highly non-linear behavior in the FEL oscillator is difficult to analyze due 
to the complexity that arises from the multiple longitudinal optical modes interacting with the 
electrons. One of the consequences of this nonlinear behavior is the trapped particle instability which 
can be analyzed through longitudinal, multi-mode analysis [Ref. 16-21]. In this chapter, a simpler 
model will be developed that limits the number of modes in the oscillator to two, which is sufficiently 
simple to study the stability of the FEL oscillator, but not excessively simple so as to be unrealistic. 
The two-mode model presented here will provide a valuable tool to study long term stability of the 


optical wave, the development of sidebands and the onset of chaos in a FEL oscillator. 


B. ELECTRON AND OPTICAL DYNAMICS 

To study the electron and optical dynamics in the oscillator with two optical modes present, 
equations vi 1.:otion for the electrons must be developed and then used to derive the driving current 
in the optical wave equation. This can be accomplished in a manner analogous to the derivation of 
the pendulum and wave equations of Chapter II, producing two-mode pendulum and wave equations. 
These two-mode wave and pendulum equations will then be used in computer simulations to explore 
stability in the FEL oscillator. 

1. Electron Dynamics and the Two-mode Pendulum Equation 

The derivation of the two-mode pendulum equation is very similar to the single mode 

pendulum equation of Chapter II, so only a sketch of the two-mode pendulum equation derivation will 


be presented. With two optical modes present, the electron equations of motion are - 


dyB) _ -2 (+E +6 x(6,+8+B)) , 
dt mc 


ay. _! aye, _ 4-1b 
dt mc Be ED : 


where E and B are the electric and magnetic fields of the primary optical field, and E, and B, are 
the electric and magnetic fields of the secondary optical field. The charge and mass of the electron 
are denoted by e and m, and ¢ is the speed of light. The undulator’s magnetic field is denoted by 


B_. Assume that the primary optical field is of the form 


E = E(cosy,-siny,0) , B = Esiny,cosy,0) , a 
and the secondary field is of the form 
E, = E,(cosy,,-siny,,0) , B, = B(siny,,cosw,,0) , 4-3 


where p = kz - wt + and py, = k, - wt + oy, k and k, are the wave numbers, w and w, are the 
frequencies, and @ and @, are the phases of their respective optical waves. Also assume the 


undulator field is helical of the form 
B,, = B,(cos(k,z),sin(k,z),0) , sa 


where k, is the wave number of the undulator’s magnetic field and z is the longitudinal coordinate 
for the long axis of the undulator. 
Substitute 4-2, 4-3 and 4-4 into the equations of motion, 4-1, and solve for the perpendicular 


motion as in the derivation of the pendulum equation in Chapter II. Use of the approximations, 
lE\(1 - B,) <1B,| and {£,1(1 - B,) <1B,| (Ref. 8] gives : 
K ‘ 
B, = 5 Coosthrsin(k2).0) 4-5 


where K = eB, /2nmc? is the undulator parameter. The energy equation is 
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4-6 
/ 


Now a relation between the electron’s phase with respect to the optical wave and 7 needs 


to be derived. As in the single mode pendulum equation, start with y~? = 1 - p? - p?. “Using the 


perpendicular motion solved for in 4-5, solve for B, yielding 
p= 1 - (i+K’) ; ; 4-7 
z 2y? 


since the electrons are relativistic. Taking the time derivative yields _ 
a - ede 
z 2 
Y Y 
Now define the electron phase with respect to the primary optical wave so that ¢ = (k +k,)z - wf. 


Taking the second time derivative of ( yields 

peat 

(k +kJc 
when solved for B,. 
Equating 4-8 and 4-9, and using the resonance condition A,/24 = y7/(1 + K?) gives 
penta | ee See 2S 4-10 

(k+kJc 2k, 2k” 7 
since kj<k. Equating 4-6 and 4-10, combining with the definition of ¢ and introducing the 
dimensionless time + = ctf/L and dimensionless optical fields fal = 4nNKeLE/y?mc? and 


jaJ = 4nNeKLEJy?mc? gives the two-mode pendulum equation 
v =f = jalcos(( +) + lalcos(kz+¥,) . ye 


As before, "°" denotes d/dt. The z dependence in the pendulum equation can be eliminated by 
using the approximation k, k, > k, and the definitionAv, = ((k-k)/k)2nN = (A, - AYA)2ZnN 


where N is the number of undulator periods. This gives 











aa 


© = |alcos({ +) + la,lcos({ + Avr +o) . ” 412 


Normally Av, is between 2x and 10x since the areas of interest are (A, - A)2nN/A * a'? which 
is the typical range of synchrotron oscillations. It is interesting to note that by defining ¢ with 
respect to the carrier wave (k,w), the symmetry of the two-mode equation is broken. ; 


2. Optical Wave Evolution and the Two-Mode Wave Equation 





The two-mode wave equation is derived in the same manner as the single. mode wave 


equation. The optical fields of 4-2 and 4-3 have the vector potential _ 


zs E. “ 
A = Z(siny,cosy.0) + —(siny,,cosy,0) 413 


Substituting 4-13 into the general wave equation 2-24 and applying the slowly varying amplitude and 
phase approximations (SVAP) to eliminate all terms with two derivatives yields the left hand side of 


the wave equation 


( = 1¢e { -2E(cosy,-sin¥.0) 


c? ot? 

2 4-14 

- -Ed(siny,cosy,0) 
c 
a. 

+ —E,(coswp,,~siny,,0) 
c 
2 

- —Fsbs(sinv cosy ,0) , 

The rapidly rotating y terms of 4-14 can be eliminated by defining four polarization vectors 
£, = (cosy, -siny,0) £, = (siny,cosy,0) as ete 


£,, = (cosy,,-siny,,0) €,. = (siny,,cosy,,0) 


to project into the wave equation producing four equations 
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Ze + E,cos(y - p,) + E,6,sin(y - y,)) = y, "é, 


2(-Bb + Esin( - ¥,) - Eyb,cos(y - ¥))=-A2J, -2, 


4-16 
2 Beosty - ¥,) - Ebsin(y -w) +B) = -27, -&, 


4x 
Cc 


2 Esincy -p,)-Edcos( ~ ¥,) - Eb) = - 


J 1° €s2 
Now the right side of the general wave equation needs to be developed to find the transverse beam 
current. Since a single particle has a transverse current, J y® -ec fi, 8% - F) (Ref. 9], giving 
J ‘ge (ecK/y)(cosk,z,sink,z,0)8(# - F,) when combined with the transverse electron motion, 4-5. 


Projecting each of the four polarization vectors into the single particle current yields 


Jy 8 = so cost hye + WBE -F) 


J, -@, = “Asingk,2 +w3e -F 

Y 4-17 
J, es, = ee cok + HbR - A) 
Jy 8s: = ech sin(kyz + WOR - F) 


Y 
Summing the single particle currents, 4-17, over all the particles, equating with the left hand side of 


the wave equation, 4-16 and multiplying each of the four resultant equations of motion with 


4nNeKL?/y’mc? yields 


|a| + {a,Jcos(Av,t +,-) - la,ldsin(Av,t +4, -6) = -Kcos(f +6)) 
lald + |d,lsin(Av,t +, - 9) + |a,|4,cos(Av,t +, -4) = ssinS +4) 
|a|cos(Av,t +o,-6) + |a|$sin(Av,t +o,-o) + la, = -j(cos(k,z + $,)) 


-|a|sin(Av,t +o,-6) + |a|dcos(Av,t +0,-) + la,16,= Jisin(k,z +y,)) 


when the dimensionless time 7 is introduced and j is the dimensionless current defined as before. 


Multiplying the second and fourth equations of 4-18 by i and then combining the first and the 
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second, and the third and the fourth equations yields the two-mode wave equations 


+ a, eft = _je-®) 


Qo 


4-19 


a,e*™ = pie Shree) 


Qo 


since k,+w, = (1 - Av,/2nN)C(t) + Av.t +,. Consistency between the two equations of 4-19 
requires Av <2xN. This eliminates the second of the equations in 4-19 so that there is only one 


wave equation to follow 


a “ a,e“"* = -jlei8) 4-20 


3. Computer Implementation of Two-mode Theory 

The two-mode wave equation, 4-20, is not immediately implementable for computer 
simulation since the two-modes need to be separated into two separate equations. To separate a froma, 
in 4-20, one simplifying approximation must be made, the gain over a single pass will be small. This 
limits the two-mode theory to FEL oscillators with low single pass gain since amplifiers have large 
gain over one pass. This approximation allows the wave equation to be integrated over time to 
eliminate one of the optical terms. Before separating the two optical terms in 4-20, a value for Av, 
must be established. Since the sideband is the mode that is modulating the carrier wave, and is 
produced by the synchrotron oscillations, it will have a frequency that is different fran the carrier 
wave by 21, the frequency of one synchrotron oscillation. Therefore Av, = 21. 

To separate the two-modes out of 4-20, integrate 4-20 over time t from 0 to I.using the 


two approximations above. This yields 


4-21 
ha = -j( ffeBrar) 


To find the change in the sideband, multiply 4-20 by e~"* and integrate over time again yielding 


ha, = [le de) an 
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Since it is easier to deal with real rather than complex numbers in computers, equations 4-21 


and 4-22 can be broken into real and imaginary parts by application of Euler’s theorem yielding 


Aa, = -j{' (costae 


Aa, = j{' (sing)de 
4-23 


Aas, = ~i{' (cos(t + 2nt))dt 


Aa, = j { Gin(G +2nt)\dt .  * 


Similarly, the two-mode pendulum equation, 4-12, can be broken into its real and imaginary parts 


= 68 +a, cost, - a,sing + a,,cos( +2nt) - a, sin(C +2nt) , 4-24 


where 5 is the tapering constant or artificial acceleration applied to the electrons after saturation in 
order to maintain resonance. Unlike an FEL amplifier, the tapering constant is always on in a FEL 
oscillator if there is to be tapering applied. . 

The program to implement the two-mode model is now straight forward. There are three 
main loops. The first, inner most, loop integrates the pendulum equation, 4-24, over all the electrons 
and keeps a sum of cos€, and cos(¢, + 2%) to be used as the averages in the wave equation. The 
next outer loop is the time loop repeated for each time step of a single pass. It integrates the wave 
equations using the averages produced by the inner most electron loop. The averages are reset for 
each time step. The outermost loop is the pass loop that controls the number of pases through the 
undulator. At the end of each pass, the optical wave is attenuated by the oscillator Q factor. 

4, Simple Oscillator Model 7 

For low-currents and weak-fields, the gain can be calculated from 2-42 as sopdcnaiinaalyG : 0.135) 

for v, = 2.6 for each pass through the undulator. For many passes, the power will beP = Pie™ 
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where n is the number of passes. The loss of optical power due to cavity losses is related to the "Q" 


factor of the oscillator so that over many passes the power in the oscillator will be 4 
P = Pyexp(0.135jn - n/Q) = Plexp(O.135jQ-In/Q . 7 4-25 | 


From this simple model, any combination of j/Q = C/ and n/Q = C2, where C/ .and ©2 are 
constants, will scale the power and phase evolutions to look identica; atter n passes. With this 


simple model, many different combinations of n, j and Q can be explored with one simulation. 


C. RESULTS OF COMPUTER SIMULATION OF TWO-MODE MODEL 
1. Comparison with Existing Simulations 

In order to establish the validity of the two-mode model and the correctness of its integration 
routine, simulations using the two-mode model over one pass are compared with existing simulations 
known to give correct results over one pass. Figure 4-1 and 4-2 compare the two-mode model with 
the single-mode model used to produce the earlier FEL simulations. In both figures, the large 
window is the phase-space evolution of the electrons with the final phase-space position plotted at 
the end of the undulator pass. Over one pass, the evolution of the second mode will not have a 
significant affect on the power, phase or electron phase-space evolution if it is started at zero power 
since it does not grow significantly. From Figure 4-1 and 4-2, it is apparent that the phase-space 
evolution is nearly identical. The optical power and phase in both simulations produce nearly the 
same final powers and phase. From this comparison it can be seen that the two-mode’s integration 
technique is probably correct. 

The two-mode model can be compared with a much more sophisticated model to 
demonstrate that the two-mode model accurately predicts the onset of a single sideband. Figure 4-3 


shows a simulation of the two-mode model for j = 5 and Q = 7 over 1000 passes. Note how 


quickly the main mode grows, saturates, and has a decaying oscillation. The sideband has a similar 
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Figure 4-1: Single-mode simulation for comparison with the two-mode model. 


pattern but uscillates about a much lower power, about 1/5 the power of the main mode. Figure 4-4 
is a multi-mode simulation where each mode is freely allowed to evolve at its natural frequency. This 
simulation has a random fluctuation with a rms value of 5a = 5 added to the initial optical field 
strength to simulate noise in the optical wave. The left-center window displays the evolution of the 
optical power, la(z,n)l, over 1000 passes through a gray scale where black is la(z)i = O and white 
represents la(z)l = 25.5 with a contour line at la(z)! = 12.8. The window has a width of two 
slippage distances with periodic boundary conditions. The slippage distance is the interaction distance 
that the electrons can couple with the optical wave. The slippage distance is the dietetics that the 
optical wave gains over the electron pulse while overtaking it during the trip down the undulator and 
is given by, NX, where N is the number of undulator periods and A is the wavelength of the main 


lobe of the optical wave. For each pass, the electrons start almost uniformly spread in phase-space 
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Figure 4-2: Two-mode model for comparison with figure 4-1. 


at each site along z, but have a small random displacement in phase where a rms spread of 
5¢ = 0.000! is added to simulate shot noise. The optical wave begins with no modulation but 
undergoes substantial gain during the first 100 passes. The bottom left window plots the gain per 
pass. By n = 50, the optical wave is already displaying some modulation and by n = 200 the 
modulation is easily visible as the fluctuations in intensity across the middle-left window. The top- 
Jeft window shows the final power at the end of the 1000th pass which is clearly modulated. The 
modulation is caused by the electrons which slip back past the optical field sites and execute 
synchrotron oscillations. The synchrotron frequency is imposed on the optical wave over many 
passes, so that only a small fraction of a synchrotron oscillation is sufficient to pass information to 


the optical wave. The bottom-right window plots the total average power in the oscillator. Note the 


power saturates at n = 100, but begins to grow slowly at n = 250 as the modulations in the main 
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Figure 4-3: Two-mode simulation with j = 5, Q = 7, for comparison with figure 4-4. 


lobe become more pronounced from the growing sideband. The sideband and the main lobe are 
clearly visible in the middle-center window where the power versus frequency is plotted with the 
same gray scale at the end of each pass. The main lobe is clearly visible as the thick, white line 
down the center of the window with the sideband approximately 2 away becoming visible at about 

= 250. The top-center window is a plot of the final power versus frequency. The ratio of power 
between the sideband and the main lobe is about 1/5 which compares favorably with the two-mode 
simulation of Figure 4-3. The right-middle window is a plot of electron phase-velocity at the end of 
each pass. It is interesting to note the phase velocities spread out as the sideband becomes more 
pronounced. The top-right window plots the final electron phase-velocity distribution. The bottom 


center window is a plot of the gain spectrum. In the multi-mode simulation, the total final power is3.2 x 10? 


which compares favorably with the total power of 3.2 x 10° in the two-mode simulation. 
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Figure 4-4: ee simulation with j = 5, Q = 7, showing the onset of a single sideband. 


Figure 4-5 and Figure 4-6 show a two-mode and multi-mode simulation similar to the 
examples above, but now Q@ = 9. The higher Q value allows the sideband to fully develop to a 
power comparable to the main lobe. In Figure 4-5, the two-mode simulation shows the sideband 
being fully developed at n = 390. Before becoming fully developed, the sideband undergoes two 
large fluctuations before reaching the power of the main lobe and stabilizing. These fluctuations will 
be examined in more detail in the next subsection. The multi-mode simulation shows the sideband 
developing much sooner, and being completely developed by n = 170 without the large fluctuations 
of the two-mode model. In the multi-mode simulation, single-mode saturation occurs at = 50, but 
growth of the total power continues as the sideband grows from n = 50 to n = 170. A second 
sideband begins to appear at n =~ 180, but never grows to a comparable power. The two-mode 


simulation has a total power of P_ = 4.4 x 10? which compares well with P = 5.9 x 10? for the multi- 
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Figure 4-5: Two-mode simulation showing sideband growth that is comparable to the main lobe. 


mode simulation. In the multi-mode simulation, the optical power is strongly modulated from the 
very beginning of the simulation because of the early growth of the sideband. 

The fast four examples were not to show that the simple two-mode model could get the 
same answer as the much more complicated multi-mode simulation, but were to show that the simple 
two-mode model is basically correct and produces results that can give insight into the development 
of sidebands. In the next subsection, the two-mode model produces results that give insight into an 
instability that leads to the onset of sidebands. Comparisons with the multi-mode simulation show 
that the instability is real and not an artificial feature of the two-mode model. It will be the simplicity 
of the two-mode model that will hopefully allow greater insight into the dynamics of the instability 


at the onset of sidebands. 
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Figure 4-6: Multi-mode simulation with early growth of the sideband to powers comparable to the 
main lobe. 
2. The Onset of Sidebands 
Figures 4-7 through 4-12 illustrate a very interesting phenomena during the development 
of a sideband as Q is increased gradually over different simulations while holding j = 5 constant. 
In Figure 4-7, where Q = 5, the sideband only grows to about 1/5 of the power of the main lobe. 
Both the main mode and the sideband saturate very early at n ~ 5Q, and the sideband has a small 
"ripple" that quickly decays away. The main lobe achieved a power of P = 2.0 x 10°. In Figure 48, 
Q = 6, and the power in the main lobe is slightly higher at P = 2.5 x 10’, but the sideband is still 
about 1/5 of the power of the main lobe. The "ripple" in the sideband decays slower than in Figure 
4-7 and the main lobe is beginning to pick up a similar "ripple". The optical phase changes rapidly. 
The "ripple" is now very pronounce in both the main lobe and the sideband. In Figure 4-9, 
Q = 7.47, the “ripple” no longer appears to decay. In Figure 4-{0, Q = 7.5, the "ripple" is a 


growing oscillation. The mean power of the sideband is still about 1/5 of the mean power of the 
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main lobe. In Figure 4-11, Q = 7.8, several interesting changes happen in the sitiddion, First, the 
sideband's average power grows and the oscillations in the sideband’s power evolution ae the main 
lobes’s power evolution twice. Second, the main lobe’s average power falls. Third, the sideband’s 
power evolution, after crossing the main lobe’s power evolution twice, oscillates to zero, then builds 
up to a steady power at nearly the same power as the main lobe. Fourth, in all the earlier 
simulations, the electrons were contained within a 27 region of phase-space so that AC $27, and 
only had to cross the periodic boundaries of phase-space once. In Figure 4-11, AC = 3n since the 
electrons are widely spread in phase-space. Finally, note the difference in the optical phase evolution. 
In the previous examples, the optical phase of the sideband and the main lobe evolved rapidly, but 
at a constant rate and together. In Figure 4-11, the optical phase evolves much slower after the 
sideband reaches steady state. In Figure 4-12, Q = 8, the sideband’s power evolution crosses the 
main lobe’s power evolution twice before it reaches a steady state. Again, the optical phase evolution 
is slower after the sideband reaches steady state, and the electrons have a large spread, 8€ < 3n. 

The oscillations and instability presented in Figures 4-7 through 4-12 are not due to 
numerical integration errors. Numerous simulations were tried with many more particles and time 
steps and the oscillations still appeared with the same shape and frequency and only slightly different 
amplitude. Near the critical Q value that causes the "ripple" to start growing instead of decaying, 
the number of particles and time steps do affect the pass that sideband develops at. The "ripple" still 
has the same frequency. The slight differences between the simulations with more time steps and 
particles are most likely due to numerical accuracy. Figures 4-13 and 4-14 show simulations with 
200 particles and 201 time steps to compare with Figures 4-3 and 4-11 which were run with 30 
particles and 30 time steps. 

The sequence of two-mode simulations in Figures 4-7 through 4-12 can be compared with 
multi-mode simulations at critical values of Q to show that certain key features exist between the two 


models. Some care must be exercised when making direct comparisons between the two-mode and 
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Figure 4-J3: A simulation of 200 particles and 201 time steps at Q = 7. 


multi-mode model since the two-mode model does not demonstrate much sensitivity to noise. Noise, 
however, has a strong influence on when the sidebands begin to develop in the multi-mode model 
and sidebands will not develop without some noise present in the multi-mode simulations. Figures 
4-4 and 4-6 of the previous subsection are far enough away from the critical values of j and Q near 
the instability, that the noise sensitivity does not appear to affect the results substantially and the 
multi-mode model accurately reproduces comparable sideband power predicted by the two-mode 
model. 

In Figure 4-15, the multi-mode simulation is run with j = 5,and Q = 7.5. The simulation 
also has a large amount of shot noise, ¢ = 0.001 which tends to cause the sideband to develop 
early. This region of j and Q is similar to Figure 4-10, so that the sideband would now have a 


growing oscillation. Since the multi-mode simulation is run to 5000 passes, there is probably time 
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Figure 4-14: A simulation with 200 particles and 201 time steps at Q = 7.8. 
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for the sideband’s power evolution to grow and overlap with the main lobe’s power oscillation as 
happens dramatically in Figure 4-11. The large shot noise probably made this interaction happen at 
an earlier time than in the two-mode simulation. This overlapping of the side and main lobe power 
evolutions could be the cause of the sideband that appears to "grow" out of the main lobe in the 
center middle window of Figure 4-15 and then "grow" back into the main lobe. This process could 
repeat itself many times as the sideband and main lobe’s oscillations slowly grow and merge as in 
the two mode model. It is also interesting to note that the electron phase velocities became more 
dispersed upon the first appearance of the sideband oscillation in Figure 4-15. This phenomenon also 
was apparent in the two-mode model. Note that the total power plot in the bottom right of Figure 
4-15 is constant after the first 1000 passes implying that the main lobe and the side lobe are 


exchanging power as the merge and re-emerge with each other. This will produce the out of phase 
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Figure 4-15: The emergence of a sideband out of the main lobe. 









oscillations that are visible in the two-mode simulation of Figure 4-10. 

Figure 4-16 is a multi-mode simulation with j = 5 and Q = 7.6 with no shot noise, but 
a large random fluctuations in the initial optical power, 5a = 5. Here, the sideband starts growing 
almost immediately which is visible by the slight grey trail that forms early just to the right of the 
main lobe in the center-middle picture of Figure 4-16. However, it is not until n ~ 280 that the 
electron phase velocity widens substantially which is when the sideband is at nearly half the 
maximum power of the main lobe. The main lobe and sideband never merge in this simulation, but 
they both have oscillations across different frequencies as they evolve. Note that the total power does 
not become constant until 1 = 400. 

The final three figures of this chapter demonstrate the simple oscillator model of section B.4 


with the two-mode model. In these examples, jQ = 39 and n/Q = 128.20513 and are constants. 
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Figure 4-16: Multi-mode simulation with the sideband developing early. 


In Figure 4-17, j = 2, n = 2500 and Q = 19.5. In Figure 4-18, j = 1, = 5000, and Q = 39. 
In Figure 4-19, j = 0.5, n = 10000, and Q = 78. Note all three simulations are similar, and the 
differences are probably due to numerical error. These figures also are similar to Figure 4-11 where, 


j = 5.x = 1000, and Q = 7.8. 


D. CONCLUSION 

Even though no quantitative theories for the development of sidebands in the FEL oscillator 
have been developed from the two-mode model, the two-mode model was instrumental in the 
discovery of the apparent instability that exists just prior to the growth of the first sideband. Together 
with the multi-mode model, the two-mode model will provide a powerful tool for the exploration of 
sideband development. It is the simplicity of the two-moue model that will probably provide the 


necessary intuitive understanding of the dynamics of sideband development to allow a simple, 
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Figure 4-17; Simple oscillator comparisons, j = 2, Q = 19.5 and n = 2500. 
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quantitative theory to be developed explaining sideband development. The two-mode model can be 
easily extended to three or more modes without excessive loss of simplicity so that studies can be 


conducted o1. the onset of chaos as well as the growth of many sidebands. 
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i 
Figure 4-18: Comparison of simple oscillator model, j = 1, Q = 39, n = 5000. 
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